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ABSTRACT 


Name of student: Paritosh Kumar Roll No: Y3 1 0303 1 

Degree for which submitted: M.Tech Department: CE 

Thesis Title: FINITE ELEMENT ANALYSIS OF MULTILAYERED REINFORCED 
GRANULAR BED SOFT SOIL SYSTEM 

Name of thesis supervisor (s): Dr. Sarvesh Chandra 

Month and Year of thesis submission: June, 2005 

The present work pertains to multi layered reinforced earth foundation system 
and uses finite element analysis to predict the settlement response of foundation as 
well as reinforcements. Strip footing and combined footing resting on multilayer 
reinforced earth are considered for load acting at the centre and at the edges 
respectively. Reinforcements are placed evenly spaced in the dense compacted 
granular bed overlying loose soil deposit and one layer of reinforcement is placed at 
the interface between the loose soil layer and dense compacted granular bed. In this 
mechanical model foundation and reinforcements are modeled as one dimensional 
beam and soil is modeled as Winkler springs. Finite element formulation for 
foundation beam, reinforcing beam and soil are presented and elemental stiffness 
matrices for the same are obtained. The matrix equation in the final form is derived 
by assembling the elemental stiffness matrices and the load vectors and solved by 
Gauss elimination method after applying boundary conditions. The model results are 
validated by comparing it with closed form solution by Maheshwari (2004) in a 
degenerated case of one layer of reinforcement. A detailed parametric study for a 
range of non-dimensional parameters is carried out and the results are obtained for 
multilayer reinforced earth foundation system. It is found that increasing number of 
layers of reinforcement beyond three does not significantly contribute to the 



settlement response of foundation beam. For three layers of reinforcements the 
reduction in maximum deflection of foundation beam are 23 % and 35 % for the load 
acting at the centre and at the edges respectively. Effects of various parameters on 
settlement responses of foundation and reinforcing beams for both loadings are found. 
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CHAPTER 1 


INTRODUCTION 


1.1 General 

With the rapid development of infrastructure in the present era there is an 
increased demand for land. Scarcity of good land and their exorbitant cost forced the 
developers to construct on the poor soils having low supporting capacity giving birth 
to the development of a new area called ground engineering. One of the ground 
improvement techniques is reinforcing the poor soil with inclusion of reinforcing 
material e.g. bamboo strip, ropes, metal strips, geosynthetics etc. to create a composite 
material with improved performance characteristics. The use of geosynthetics for soil 
reinforcement is very extensive and increasing annually. The simplicity of the basic 
principles of geosynthetic reinforcements and the economic factors over the 
conventional approach makes it very attractive to the civil engineers. A large number 
of experimental and analytical studies have been carried out in the last three decades 
in the area of footings resting on soil reinforced with geosynthetic reinforcement. The 
analysis of the interaction between stmctural foundations and supporting soil media is 
of fundamental importance to geotechnical engineering. Results of such analyses 
provide information, which can be used in the design of foundations and in the 
analysis of stresses and deformations within the supporting soil medium. Various 
analytical, numerical and experimental methods are used for the analysis of 
foundations on reinforced ground. This present study analyzes the foundation beams 
resting on earth reinforced with multilayer geosynthetic reinforcement by finite 
element method with lumped parameter approach for the settlement response of 


foundation beam and reinforcements. 



1.2 Scope and organization of the present work 

A brief review of theoretical works on the . behavior of foundations on 
reinforced soil is presented with conclusive remarks in chapter 2. 

From the review of the literature, it has been observed that a very few works 
have been done on reinforced earth foundation system taking into consideration the 
bending stiffness of reinforcement either for single layer reinforcement or for multi- 
layer reinforcement. Thus it is felt that there is a need to undertake such a study and 
the same is presented in the following chapters. 

Problem formulation and method of analysis is presented in chapter 3. In this 
chapter statement of the problem is presented. The formulation for foundation beam, 
reinforcing beam, soil and the effect of frictional force and prestressing force in beam 
in finite element form is presented subsequently. 

Results and discussions of the present study are presented in chapter 4. Effects 
of various parameters on the settlement behavior of the foundation and reinforcement 
are discussed. 

Chapter 5 deals with the conclusions of the present study and finally future 
scope of the present work is presented. 



CHAPTER 2 


LITERATURE REVIEW 


2.1 General 

Ground improvement techniques are one of the most interesting and widely 
studied topics in geotechnical engineering. Reinforcement of soft soils to improve the 
bearing capacity and reduce settlements of structures built on them is a very efficient 
method and is becoming increasingly common. Geosynthetic reinforcements have 
been proved to be very efficient and cost effective in ground reinforcement. There has 
been a need to study the geosynthetic reinforced soil foundations to predict the 
behavior of foundation and soil. Several studies have been done on the geosynthetic- 
reinforced earth in the past three decades to find the effective use of such foundation 
systems. The brief review given in the following section refers to the literature related 
to the theoretical works carried out in the area of reinforced earth to study the 
response of the foundation on reinforced soils. 

2.2 Geosynthetic reinforced soil foundation system 

Binquet & Lee (1975) analyzed the bearing capacity problem of a surface strip 
footing on a granular soil reinforced with horizontal layers of tensile reinforcement. 
The failure surface was determined as the locus of the points of maximum shear stress 
at different depth based on the elastic theory of stresses i.e. by applying the 
Boussinesq equation. It was arbitrarily assumed that the tie force per layer varies 
inversely with the number of layers, N in the foundation due to lack of definitive data. 
Idealized approximate cost analysis indicated that if corrosion was neglected, savings 
up to 100% can be considered over the cost of conventional foundation. As 



Boussinesq equation was found to be unsatisfactory for determining stresses in the 
case of the foundation on granular soil overlying a finite sized pocket of soft material, 
linear elastic finite element program was used for the analysis. 

Brown & Poulos (1981) presented an analytical model to investigate the 
increase in bearing capacity and stiffness of a foundation due to the placement of 
reinforcement in the homogeneous soil layer, and soil bed overlying a cavity. The 
analysis considered various components of the reinforced soil separately and 
incorporated an elasto-plastic Mohr-Coulomb failure criterion, equivalent-width 
elastic reinforcement strips with zero flexural rigidity transmitting axial forces only, 
and a model for reinforcement-soil bonding which allowed for slippage governed by 
Mohr-Coulomb failure criterion. It was concluded that the improvement of the 
performance characteristics of the reinforced soil bed depended both on the number of 
reinforcing layers and also on the concentration (surface area per unit width of the 
footing) of the reinforcement, the latter being the major dominant factor as it was 
observed that the limit state of the reinforcement-soil bond was reached at an early 
stage. However, the study was only conducted for four layers of reinforcement due to 
lack of resources and thus needed further investigation about the effect of number of 
reinforcing layers on the performance characteristics of the reinforced soil bed. 

Andrawes et al. (1982) analyzed soil-geotextile systems by using finite 
element method. They described the nature of the elements used to represent the soil- 
geotextile systems for the purpose of predicting the stress-strain behavior. Studies 
were done with a footing resting on dense sand without or with geotextile placed at 
various depths. It was observed that if the soil properties could be correctly 
represented in the soil elements, good correlations up to 85% were obtained between 
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predicted and measured data. However, it was suggested that the finite element 
method was inappropriate for the soils in which local failures occur. 

Madhav & Poorooshasb (1988) proposed a new foundation model element - 
the rough membrane to represent the response of the geofabric. Combining this 
element with the Winkler springs and Pasternak shear layers to model respectively the 
soft soil and granular fills, the foundation model was proposed for geosynthetic- 
granular fill-soft soil system. Analysis of results at small displacement indicated the 
effect of granular fill to be more significant than that of the membrane. It was also 
observed that the effect of membrane increased with the load or decreasing soil 
stiffness. 

Poran et al. (1989) presented a design procedure based on finite element 
analysis, which included a visco-plastic model for soils and special visco-elastic 
membrane elements to model geogrid behavior. It was suggested that the footing 
settlement is the governing design criteria, and therefore granular layer thickness and 
footing size should be determined based on allowable settlement, especially with 
geogrid reinforcement. The proposed design procedure for continuous and 
axisymmetrical footing was especially beneficial where small settlements were 
encountered. 

Ghosh & Madhav (1994a) developed a simple mathematical model to account 
for the membrane effect of a reinforcement layer on the load settlement response of a 
granular fill-soft soil foundation system. The nonlinear loading pressure-settlement 
response of the granular fill and soft soil, respectively, for the plane strain condition 
were incorporated in the formulation. Parametric studies for a uniformly loaded strip 
footing showed that the membrane action of reinforcement caused reduction in the 
settlement below the footing which is over and above the settlement due to granular 



fill. The type of foundation considered consisted of a layer of frictional fill over soft 
clay. In the reinforced case, the reinforcement was placed within the granular fill and 
at the fill-clay interface. The normal stresses at the reinforcement interfaces were not 
much affected by their presence and the internal stability of the reinforced soil system 
might be impaired by the tensile rupture of the reinforcement or by lack of adherence 
of the soil at the reinforcement interface. However, in the reinforced case, analysis 
with multi-layer of reinforcement was not undertaken. It was observed that the 
settlement behavior significantly improved with respect to the stiffness of granular 
fill, when the soil was softer, and with respect to the interfacial friction, when the fill 
material was stiff. 

Ghosh and Madhav (1994 b) developed a mathematical model for the analysis 
of a reinforced foundation bed by incorporating the confinement effect of a single 
layer of reinforcement. It was quantified in terms of the average increase in confining 
pressure due to the reinforcement from which modified shear stiffnesses of the 
granular soil surrounding the reinforcement were obtained. The parametric studies 
were sensitive to the methods used for obtaining the modified shear stiffness of 
granular fill. The confinement effect was more pronounced when the shear stiffness of 
the granular fill was large. The modified shear stiffness below the centre of the 
footing increased by two to five times the initial values of shear stiffness. 

Ghosh and Madhav (1994 c) proposed a new model for a reinforced shallow 
foundation bed by incorporating the ‘rough membrane’ element for single layer 
reinforcement. Mechanics of the rough membrane element with the assumption of 
horizontal shear stress transfer at the soil/reinforcement interfaces were explained and 
formulated. The model was generalized by incorporating nonlinear response of the 
soft soil and of the granular fill under plane strain loading conditions. The membrane 


effect was quantified in terms of the improvement in the load-settlement responses of 
the composite foundation system over the unreinforced one. Parametric studies 
indicated that reinforcement while in tension spreads the load over a larger area, 
leading to a reduction in the settlement beneath the footing. 

Shukla & Chandra (1994a) studied the effect of prestressing the geosynthetic 
reinforcement on the settlement behavior of geosynthetic-reinforced granular fill-soft 
soil system. The foundation model element had a prestressed rough membrane 
embedded in a granular layer. Parametric studies revealed that the settlement 
reductions within the loaded regions were observed even for low prestress in the 
reinforcement. It was observed that the improvement in settlement response increased 
with the increase in prestress in the geosynthetic reinforcement within the loaded 
region and was most significant at the centre of the loaded footing, which reduced 
differential 'settlement. On soft soils or at higher load intensities more prestress was 
required in the geosynthetic reinforcement to achieve a particular settlement reduction 
compared with stiff soils or lower load intensities. The influence of Prestressing 
increased with the width of the reinforced zone up to UB ratio of 2.0. The effect of 
prestressing was significant at lower values of shear parameter of granular fill and 
interface friction coefficients at top and bottom of geosynthetic reinforcement. It was 
also observed that the mobilized tensile force in the geosynthetic layer increased at all 
locations with the increase in prestress but the increase was maximum at the edges of 
the footing. It was concluded that prestressing the geosynthetic reinforcement was a 
significant ground improvement technique to enhance the settlement characteristics of 
the soft soils where the membrane effect of reinforcement is required. 

Shukla & Chandra (1994b) proposed a foundation model to incorporate 
compressibility of the granular fill by attaching a layer of Winkler springs to the 
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Pasternak shear layer. The effect of compaction of granular fill in terms of lateral 
stress ratio was also considered. The parametric studies indicated that the 
consideration of the compressibility of the granular fill resulted in a significant 
increase in the settlement of the reinforced foundation soil. The model used was fairly 
general in nature as it accounted for the compressibility of the granular layer, in 
addition to the effects of prestressing of geosynthetic reinforcement and 
overconsolidation (through compaction) of the granular fill. The model was capable 
of incorporating the effect of any degree of compressibility of the granular fill by 
proper selection of its modulus of subgrade reaction. Finite Difference scheme was 
employed to produce the numerical solution. Non-dimensionalized results were 
produced for practical applications. The effect of the compressibility of the granular 
fill on the load settlement characteristics of geosynthetic-reinforced granular fill-soft 
soil system was observed to be significant at all the load intensities except for very 
small values of the load. With the increase in the compressibility of the granular fill, 
the order of increase in the settlement was more at the edges than in the centre of the 
loaded region. For any degree of compressibility of granular fill, the geosynthetic- 
reinforced granular fill-soft soil system behaved as much as a stiffer system at higher 
load intensities. The order of increase in the settlement with increase in the 
compressibility of the granular fill at any location within the loaded region was 
greater for lower values of shear parameter. 

Shukla & Chandra (1994c) described a mechanical model for idealizing the 
settlement response of a geosynthetic-reinforced compressible fill-soft soil system, by 
representing each subsystem by commonly used mechanical elements such as 
stretched, rough, elastic membrane, Pasternak shear layer, Winkler springs and 
dashpot. The model considered simultaneously several factors governing its behavior 



9 


such as the compressibility of the granular fill, the time-dependent behavior of soft 
soil and prestress in the geosynthetic reinforcement, besides the material and soil- 
geosynthetic interface characteristics. The response function of the model had been 
derived for strip loading in plane strain conditions. The resulting equations, in 
nondimensional form were solved by an iterative finite difference method. It was 
observed that the compressibility of the granular fill had an appreciable influence on 
the settlement response of the geosynthetic-reinforced granular fill-soft sol system as 
long as the stiffness of the granular fill was less than approximately 50 times that of 
soft soil. It was also concluded that the pre-stressing of the geosynthetic 
reinforcement and the compaction of the granular fill resulted in the settlement 
reduction of reinforced fill-soft soil system. 

Gharpure (1995) studied the finite element analysis of geotextile reinforced 
sand bed subjected to strip loading. The soil was modeled as plain strain element and 
reinforcement as one dimensional line element with zero thickness. The confinement 
and string effect of reinforcement was considered along with the slip between 
geotextile and surrounding soil and gave the optimum depth of placement for 
geotextile. 

Yin (1997) proposed a new one-dimensional mathematical model for 
modeling geosynthetic-reinforced granular fills over soft soils subjected to a vertical 
surcharge load. The geosynthetic reinforcement consisted of a membrane (geogrid, or 
geotextile) placed horizontally in engineered granular fill, which was constructed over 
a soft soil. The proposed model was mainly based on the assumption of a Pasternak 
shear layer. A new approach consisted of incorporating a deformation compatibility 
condition into the model. The compatibility condition eliminated the need for two 
uncertain model parameters that were required to calculate the shear stresses between 
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the top and bottom granular fill layers and the geosynthetic layer. The compatibility 
condition also provided the possibility to include the geosynthetic stiffness in the 
model. Comparisons were made with a two-dimensional finite element model and 
three other one-dimensional models. The proposed model could also consider the 
nonlinear and elastic-plastic behavior of the granular fill. 

Yin (2000) performed a study of reinforced beam on elastic foundation. 
Governing ordinary differential equations were derived for a reinforced Timoshenko 
beam on an elastic foundation. An analytical, solution was obtained for a point load on 
an infinite Timoshenko beam (TB) on elastic foundation. Special attention was drawn 
to the location, tension, and shear stiffness of reinforcement and its influence on 
settlement/deflection of the beam and reinforcement tension force. The TB model was 
suggested to be particularly useful in modeling a reinforced beam with or without 
considering the reinforcement shear stiffness. Practical applications of the TB model 
were recommended in modeling geosynthetics/fiber-glass reinforcement of 
foundation or pavement. A finite element (FE) model was established for the same 
infinite beam problem. Comparisons were made with ±e FE model and PB model 
(Winkler model based on pure bending theory). It was observed that the TB model 
generally predicted larger settlement and less tension force in the distance at or near 
the point load, as compared to the PB model. However, comparable results were 
obtained with the FE model. It was suggested that the TB model was superior to the 
Winkler model in considering the shear deformation of the beam and the 
reinforcement. The TB model could consider the shear stiffness of reinforcement, 
especially for the geosynthetics, which have almost zero shear stiffness due to very 
thin thickness. However, the TB model suffered from the same limitations as that of 



the Winkler model. The TB model also suffered from the limitation due to the 


assumption of zero vertical strain on the beam. 

Fakher and Jones (2001) presented a numerical simulation to model a layer of 
sand overlying a layer of geosynthetic reinforcement and super soft clay. It was 
suggested that bending stiffness should be considered while having earthworks on 
super soft clay. It was suggested that the structural behavior of the reinforcement over 
super soft clay was like a stiff plate supported by the vertical reaction force from the 
clay. If the plate is anchored horizontally, it will present a high resistance to the 
vertical load. The effect of the reinforcement pullout resistance was to provide a 
horizontal restraint for the plate. It was concluded that the higher the reinforcement 
bending stiffness, the higher the bearing capacity of the system. It was observed that 
the reinforcement with high bending stiffness reduced the surface heave, which was 
distributed over a wider area. With very stiff reinforcement, the maximum heave did 
not occur at the surface. The relative importance of the bending stiffness of the 
reinforcement was observed to be in particular cases. When the ratio of depth of 
placement of reinforcement to the footing width is small, the increase of the bearing 
capacity ratio due to the bending stiffness is significant. When the footing width 
increases, the effect of bending stiffness decreased sharply. It was concluded that the 
bending stiffness was not important unless the underlying clay layer was not in a 
super soft state. 

Maheshwari (2004) analyzed a model to predict the flexural response of beam 
resting on dense compacted soil overlying loose soil deposit with one layer of 
geosynthetic reinforcement such as geogrids at the interface of soil layers. Both the 
foundation and reinforcement were idealized as beams with smooth surface 
characteristics. The lower poor strata and upper dense soil were modeled using 
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Winkler springs of different stiffness. Effect of the depth of placement of 
reinforcement was considered by taking surcharge on reinforcement. The governing 
differential equations for the response of the beams were derived and closed-form 
were obtained subjected to appropriate boundary and continuity conditions. The 
model results showed very good agreement with the solution provided by Hetenyi for 
infinite beams on elastic foundation in a particular case. Practically no change was 
observed in the normalized deflection of the upper and lower beams when the 
normalized length ratio of beams exceeded 1.5 for the range of parameters considered. 
The normalized depth of placement of the lower beam had a significant effect on the 
deflection profile of both foundation and reinforcing beams. The maximum 
normalized deflection of the upper and lower beam increased by 37% and 45% 
respectively as the normalized depth increases from 0.5 to 1.5 for load acting at the 
centre of foundation beam. For loads acting at the edges of foundation beam the 
maximum normalized deflection of the upper and lower beam increased, by 14-15 % 
and 25% respectively as normalized depth increases from 0.5 to 1.5. The relative 
flexural rigidity of the beam, R affected the deflection at the edges of the beam more 
than at the centre. The maximum normalized deflection of the upper beam decreased 
by 11% as R increases from 50 to 1 while for lower beams it decreased 55% for a 
reduction of R from 50 to 1. The relative stiffness of soils, r, had a significant 
influence on the normalized deflection of upper and lower beam. The maximum 
normalized deflection decreases by 84 - 85 % and 32 - 33 % for upper and lower 
beam respectively as r decreases from 20 to l.The load acting on foundation beam 
had a significant influence on the deflection of beams. The normalized deflection of 
the beams could be reduced to the extent of 70-75% at the centre and by more than 
95% at the edges of the beams. 
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2.3 Motivation of present work 

From the review of literature, it is noticed that a few studies have been carried 
out on reinforced earth foundation system taking bending stiffness of reinforcement in 
to consideration. First Fakher and Jones (2001) considered the bending stiffness of 
geosynthetic reinforcement and concluded that bending stiffness of reinforcement has 
significant influence on the bearing capacity of the reinforced earth. Maheshwari 
(2004) have also considered the bending stiffness of reinforcement and analyzed the 
foundation beam on reinforced soil with one layer of reinforcement and gave the close 
form solution for the settlement response of foundation and reinforcement. The 
present study is an extension of the work done by Maheshwari (2004) for multilayer 
reinforcement and FEM is chosen as method of analysis to predict the settlement 
response of foundation as well as reinforcements. 



CHAPTERS 


PROBLEM FORMULATION AND METHOD OF ANALYSIS 

3.1 General 

Finite Element Method (FEM) is well recognized as a very powerful and 
versatile numerical tool for solving stress-strain and displacement related 
problems. FEM is specially suited for continuum problems like modeling of soil, 
but it is also well suited and accurate with lumped parameter approach, if 
parameters and size of elements are properly chosen and are well organized. In 
this thesis lumped parameter approach is followed. From the review of literature 
in chapter 2 it is observed that reinforcements like geogrids and geocells may 
have bending stiffness also. Maheshwari (2004) modeled foundation and 
reinforcement as beam and considered one layer of reinforcement at the 
interface between the dense compacted granular bed and loose soil deposit and 
gave close form solution for the settlement response of foundation as well as 
reinforcement. In the present study the work done by Maheshwari (2004) is 
extended for multilayer reinforcement and FEM is chosen as method of analysis 
to predict the settlement response of foundation and reinforcements. Statement of 
the problem is presented first then in the following sections, the finite element 
formulations for foundation beam, reinforcing beam and soil are presented and 
elemental stiffness matrices are formed for the same. Friction and prestress in 
reinforcement beam are treated as axial force in beam and the formulation is 
given for the same. Then the finite element mesh details and connectivity of the 
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various components are discussed. Finally assemblage of elemental stiffness 
matrices and load vectors to derive finite element equation in matrix form, 
application of boundary conditions and solution of FEM equation are described. 

3.2 Statement of the problem 

The cases of strip footing of length 21i acted up on by a concentrated load 
Qt at the centre is shown in Fig. 3.1 (a) and a combined footing of length 21i acted 
upon by concentrated load Qt at the edges is shown in Fig. 3.1 (b). Footings are 
idealized as elastic beams with flexural rigidity Eili resting on dense compacted 
soil of depth h placed over loose soil deposit. Geosynthetic reinforcements 
(geogrids or geocells) are idealized as beams of length Ih with flexural rigidity 
E 2 I 2 and are placed at the interface of soil layers as well as in dense compacted 
soil layer with equal spacing h/3 as shown in Fig. 3.1 (a) and 3.1 (b) for strip 
footing and combined footing respectively (Figs, are drawn for three layers of 
reinforcements for simplicity). The strong soil layer and underlying poor soil are 
idealized as Winkler springs of stiffness ki and fo. The unit weights of upper and 
lower layer of soil are yi and y 2 respectively. Fig. 3.2 shows the idealized model of 
above described reinforced earth foundation systems, which represents strip 
footing by considering the load Qt at the centre, and combined footing can be 
represented if the loads Qt are considered at the edges. The prestess force in 
foundation beam is 7} and prestress in reinforcing beams are Ti, T 2 , and 7j for 
layer - 1 layer - 2 and layer - 3 respectively and these forces are zero in 
absence of prestress forces. In this present study foundation beam is assumed to 
be smooth and frictional forces of tensile nature are developed in the reinforcing 
beam resulting from interfacial friction between reinforcing beam and 
surrounding (shown in Fig. 3.2). Loading on foundation beam is assumed to be 
symmetric and valid for strip footing and combined footing. 
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Q. 

i 



/2/3 Reinforcing beam, layer -2 (E2I2) Dense soil (ki,Yi) 



Reinforcing beam, layer - 3 (E2I2) Poor soil (k:2,Y2) 



Hard stratum 



Fig. 3.1: DeHnition sketch of the problem: (a) strip footing, (b) combined footing 
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Fig. 3.2: Idealization of the problem 


3.1 Formulation of foundation and reinforcing beams 


Foundation and reinforcement beams are modeled as a typical two noded one 


dimensional beam element. Displacement u and rotation — are considered as degrees 

dx 


of freedom at each node as shown in Fig. 3.3. Stiffness matrix for such beam element 


is well established and is used in the present study. In this section formulation for a 


typical one dimensional beam element (as shown in Fig.3.3) is presented. All 


equations in this section are taken from Tirupathin et al. (2003). 
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e 


^ du'^ 
dx 

J\ 


- 1=^1 

nodes 


u‘ 


du 

ydx j 


1=4 


; Local node 

: Elemental degrees of freedom 
: Coordinate of local node 
: Natural coordinates of local 


Fig. 3.3: A typical beam element 


3.3.1 Elemental beam stiffness matrix 

Elemental stiffness matrix for the beam element shown in Fig. 3.3 is given 


as 


[^r=f;[£/j5}'{5rk 

X\ 


(3.1) 


where 


[k] 


is the elemental stiffness matrix for the considered element ‘e’ 


w = 


dx^ 


(3.2) 


is the row vector contains the derivatives of Hermitian shape functions 
in local X coordinate. 

E is Young’s modulus of elasticity of the considered element. 


/ 


iz is moment of area of the considered element. 


{^fis 


the transpose of 


{5}’. 
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3.3.2 Numerical Integration for obtaining elemental stiffness matrix 

To numerically integrate the equation (3.1), equation is first converted to 
natural coordinates and then it is integrated by Gauss quadrature numerical 
integration scheme as given below. 

Equation (3.1) is converted in natural coordinate as 




8 £/ 


ZZ 






Gauss quadrature numerical integration scheme is 






1 / 






m=l 


ZZ 








(3.3) 


(3.4) 


where, 

is number of gauss points used. 

Gauss points used are -0.577350269189626, 0.577350269189626. 


is the weight corresponding to gauss points are 1.0, 1.0. 


E’ = t E-Ef 


2 



(3.5) 

(3.6) 


(3.7) 


1=1 

where L] are Lagrangian shape function in natural coordinates for nodes i = 1, 2 


for typical element ‘e’ is given as 
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UU - ^ j) 

^ (3.8) 

n ^ ) 

j = 1 


Shape function derivatives are given as {s}* = 


‘‘'{nY 

di' 


where, 

Hermitian shape functions in natural coordinates {n}' are given as: 


= 

■''' 2 1-1 


dL, 




LI 




(3.9) 


(3.10) 


(3.11) 


where 


L, are Lagrangian shape functions for nodes i = 1, 2 of element e. 
le is length of beam element. 


Shape function derivatives S‘ 


d!ML 


in natural coordinate scheme are 


i;'=|(-2+6f) 

S \ =^(2 + 6 fl 
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3.4 Formulation of soil 

Soil is modeled by Winkler springs and is found to be very effective in 
lumped parameter approach. In this section elemental stiffness matrix for soil 
element is found. A typical two noded line element is chosen for Winkler spring 
as shown in Fig. 3.4. Displacements at nodes («) are considered as degree of 
freedom for the spring element. Since each node have two degrees of freedom 

displacement («) and rotationf— Ifor the foundation and reinforcing beams, in 



A typical elemental stiffness matrix for spring element is given as 



0 -k 0 
0 0 0 
0 it 0 
0 0 0 


where, k represents the spring stiffness of corresponding soil. 
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3.5 Formulation for the effect of prestress force and friction force 

In this thesis, prestress and frictional forces are considered as axial force 
in foundation and reinforcement. As foundation and reinforcement are modeled 
as beam, prestress and friction effect are modeled as beam with axial force. All 
the equations given in this section are taken from Clough and Penzien (1993). 
When a beam element is subjected to an axial force, the modification is made to 
stiffness coefficient Kij by geometric stiffness kaj, which is defined as the force 
corresponding to nodal coordinate i due to a unit displacement at coordinate j 
and resulting from the axial forces in the structure. These coefficients are 
evaluated by application of principle of virtual work. Consider a beam as shown 
in fig. 3.3 but now subjected to a distributed axial force per unit of length Pa(x), 
as shown in Fig. 3.5 (a). In the sketch in Fig. 3.5 (b), the beam segment is 


subjected to a unit rotation at the left end, 


^dXj, 


=1. 



X 



Fig. 3.5: (a) variation of axial force (b) axially loaded beam with real rotation 
and virtual translation of node (taken from Clough and Penzien (1993)) 
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By definition, the nodal forces due to this displacement are the 
corresponding geometric stiffness coefficients, for example kon is the vertical 
force at left end. If now, this left end is subjected to a unit displacement m'= 1, 
the resulting external work is 

We=kGI2Ul (3.12) 

We=kG12 

Since u" = 1 


The internal work done during this virtual displacement is found by 
considering a differential element of length dx from the beam in Fig. 3.5 (b) and 
shown in enlarged form as in Fig. 3.6. 

Pa(x) 



Fig. 3.6: Enlarged elemental length dx of beam. 

By similar triangles, it may be seen from Fig.3.6. 

d{Se) _ du 
d[5u)~~dx 

d{Se) = ^d{&i) 
d{5e) = —5{du) 



(3.14) 
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The work done by the axial force Pa(x) during virtual displacement is 
dW,=Pa(x)de 


(3.15) 


Where 5e represents the relative displacement experienced by the normal force 
Pa(x) acting on the differential element during virtual displacement. 


dW, ^ P{x)—d—]dx 
' ^ ' dx \ dx 


(3.16) 


The displacements are expressed in terms of shape functions and integrating 
finally gives 




(3.17) 




dN^{x)dN,{^ 
dx dx 


(3.18) 


Comparing internal work done with external work done this geometric 
coefficient is found to be 


‘e 


dN2{x) dN^{x) 

( 

dx dx 


(3.19) 


Or in general 


I p.U) 


dN j{x) 


(3.20) 


For computer application, above integral is converted to natural coordinate then 
numerical integration is carried through Gauss quaderature numerical 
integration scheme. 

In natural coordinate system above integral equation (2.24) is 


(3.21) 



Numerical integration is carried through gauss quadrature as 



The distributed load on foundation beam or the surcharge load on 
reinforcing beam is considered as distributed load over the beam element as 
shown in Fig. 3.7 and is converted to nodal forces and moments with the help of 
Hermitian shape functions {A/'}* as shown below. 



(b) 


Fig. 3.7: (a) A typical beam element loaded with UDL and (b) their equivalent 
loads and moments (taken from Paz (1987)) 

The elemental load vector for beam element with distributed load is given by 

{/? = -yj wix){N jtir 

^ 0 


(3.24) 
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In natural coordinate above equation becomes 

{/}' (3.25) 

^ -1 

Above equation is numerically integrated by Gauss quaterature numerical integration 
scheme. Global load vector {F} is obtained by assembling the elemental load vectors 
for foundation beams and reinforcing beams then adding the concentrated load vector 
[F*] to the assembled load vector (given in the section 3.8.2). 

[F*] is given by 



(3.26) 


where 

W 2 i.i are the concentrated forces. 

M 2 i are the concentrated moments. 

i = 1,2,3 - • -are global node number at which concentrated load acts. 

3.7 Mesh details and connectivity of various components 

In this section finite mesh details i.e. global node numbering and element 
numbering for beams (representing foundation beam and reinforcing beam), springs 
in strong soil and for bottom spring representing poor soil are presented, basic 
concepts are referred from Cook (1981). Then connectivity matrices are formed for 
each component i.e. beam elements and spring elements and are used in assemblage 
of components. Finite element mesh details used to form connectivity matrices is 
shown in Fig. 3.8. 



27 



Y 


Fig. 3.8: Finite element mesh details 

1 , 2 , 3 ... are Global node numbers. 

Q] [ 2 ] [U are top spring element numbers. 

[U S Si are beam element numbers. 

[D El m bottom spring element numbers. 

Connectivity matrices for various components are given as: 

'1 5 

5 9 

Beam connectivity matrix [c * ]= 9 13 

2 6 


1 2 
5 6 

Top springs connectivity matrix [c^ ' ] = 9 10 

13 14 


Bottom spring connectivity matrix 
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3.8 Assembly of elemental matrices 


After getting the elemental stiffness matrices, elemental load vector for 
distributed load and concentrated load vector, it is required to assemble all the 
elemental stiffness matrices and the load vectors in order to get the final FEM 
equation. This section focuses on the assembling of various elemental matrices. 

3.8.1 Assembly of elemental stiffness matrices 

Elemental Beam stiffness matrix and elemental spring stiffness matrix are first 
expanded to the size of global stiffness matrix of size 2nx2n (n is the total number of 
global nodes) using simplified assembly relations given below. 




when r = Cep and s = Ceq 


=K2j 


(3.27) 


= 0 otherwise 


where, 

\_KY represents expanded elemental matrix elements for beam and spring elements. 


w represents unexpanded elemental matrix elements for beam and spring 
elements. 

Cep, Ceq are connectivity coefficients of corresponding connectivity matrix. 
p, q are corresponding local node numbers to beam and spring elements. 

/•, are corresponding position in global stiffness matrix. 


e represent the element considered. 
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Global stiffness matrix is obtained by adding all the elemental stiffness 
matrices as in equation (3.28). 




e=\ 


e=l 


e=l 


e=l 


(3.28) 


where, 

[K] is the global stiffness matrix. 

]* is expanded elemental stiffness matrix of foundation beam 

[ATJ* is expanded elemental stiffness matrix of reinforced beam 
[Kk, ]‘ is expanded elemental stiffness matrix of top soil sub grade modulus 
is expanded elemental stiffness matrix of lower soil sub grade modulus, 
w/ is number of elements of foundation beam. 


yig is number of elements of reinforcement beams. 

k 

fig' is number of elements of springs representing compact dense soil. 

flg^ is number of elements of springs representing poor soil. 

nl is the number of layers of reinforcements. 

3.8.2 Assembly of elemental load vectors 

Elemental load vectors are first expanded to the size of global load vector 
of size 2n using simplified assembly relations as given below. 

when r = Cep (3.29) 


{F2.J ={f 2 , J 
{F'2rY ={f2,} 


= 0 


otherwise 
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where, 



represents expanded load vector elements for beam elements. 


{/} represents elemental unexpanded load vector for beam elements. 

Cep, are connectivity coefficients of corresponding connectivity matrix. 
p is corresponding local node numbers to beam and spring elements. 
r is the position of element in global matrix. 
e represent the element considered. 

Global load vector is found by adding all the elemental load vectors and the 
concentrated load vector as shown below. 


e=\ e=\ 

where, 

{F}is the global load vector. 

is expanded elemental load vector of foundation beam. 

y is expanded elemental load vector of reinforcement beam. 
{F*} is the concentrated load vector. 

/ is number of elements of foundation beam, 
is number of elements of reinforcing beams. 

3.1 Final FEM equation 

FEM equation in matrix form is given by 


(3.30) 


(3.31) 


where 
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[isT] is the global stiffness matrix. 

(3.32) 

where 

Ui, are the displacements of beam at global node i=l,2,3-.-n. 

(du^ 

— are the rotations of beam at global node i=l,2,3...n. 

[dxj,. 

n is the total number of global node. 

3.10 Application of boundary condition 

Since the geometry and loading of the model considered is symmetrical, 
the rotation of foundation and reinforcing beams at the centre is zero 
corresponding elements in global load vector is made zero and corresponding 
elements in global stiffness matrix are made 1. 

3.11 Solving the FEM equation 

After applying boundary conditions the FEM equation is solved by gauss 
elimination method. 

Above formulations and analysis described in this chapter are used to 
analyze the problem considered in the present study and results from the 
analysis are discussed in chapter 4. 


{w} is the global displacement vector and is given by 


du^ 


^du^ 



CHAPTER 4 


RESULTS AND DISCUSSIONS 


4.1 General 

Following the statement of the problem and problem formulation by FEM as 
described in chapter 3, a computer program in language C is developed. The validity 
of the program and the convergence of the solution are checked first and then detailed 
results for parametric study are obtained. All the computations have been carried out 
on a personal computer with 2400+ AMD processor and 256 MB RAM. Since the 
problem considered is symmetrical with respect to geometry as well as in loading 
condition, only half portion of the system is studied and presented. 

4.2 Convergence study and Validation of the model 

The problem consists of foundation beam and the reinforcing beam having 
equal length with one layer of reinforcement subjected a concentrated load acting at 
the center is solved for three discretizations of 40, 80, 160 beam elements (for each 
foundation beam and reinforcing beam). The parameters used are /;=16m, Zr= 1, R = 
10, r = 5, Rj= 1, ki = 50 MN/m^,y = 0 and 2= 100000 N. For all three discretizations 
the normalized deflections at the centre and at the edge of foundation beam as well as 
reinforcing beam are tabulated in table 4.1. It is observed that the values of 
normalized deflection at centre and at edge almost converge for 160 elements. 
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Table 4.1: Values of deflection for three discretizations 



Normalized deflection 

Normalized deflection 


at the centre 

at the edge 

Discretization 

Foundation 

Reinforcing 

Foundation 

Reinforcing 


beam 

beam 

beam 

beam 


(yi ) 

Cva ) 

(yj ) 

Cvz ) 

40 elements 

1.22551 

0.0999994 

2.30476x10'^ 

1.94588x10'^ 

80 elements 

1.27784 

0.104006 

2.60460x10'^ 

2.19902x10'^ 

160 elements 

1.30571 

0.106116 

2.77631x10'^ 

2.34401x10'^ 


The present FEM model results are validated by comparing it with the closed 
form solution obtained by Maheshawari (2004) for both foundation beam and the 
reinforcing beam having equal length with one layer of reinforcement subjected to 
concentrated load acting at the center. While validating the program three hundred 
beam elements are taken for each foundation and reinforcing beam. The results from 
the present study for typical values of parameters are found to be in exact agreement 
with the results obtained by Maheshwari (2004) as shown in Fig 4.1. While validating 

. y E I 

the model results, the deflections of reinforcing beam are normalized as yj = ^ ^ 3 ^ 

Qt^2 


(same as Maheshwari (2004)). 
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- 4 - - Deflection of foundation beam in 
present study 

— o — Deflecton of reinforcing beam in 
present study 

- - X* - ' Deflection of foundation beeam by 
Maheshwari (2004). 

— I — Deflection ofreinforcing beam by 
Maheshwari (2004). 


Normalized distance from centre (z) 

Fig. 4.1: Comparison of deflection proflle of beams with 


Maheshwari (2004) 


4.3 Parametric study 

The model after validation for one layer of reinforcement has been used to 
study the effect of multilayer reinforcement. The effects of various parameters and the 
effect of number of layers of reinforcement on settlement response of foundation 
beam and reinforcing beams are presented in the following sections. Maheshwari 
(2004) concluded that there is no significant influence on deflection of foundation 
beam when Zr is more than 1.5 for one layer of reinforcement and the same was 
observed for multilayer reinforcement hence in this present study for all parametric 
values the ratio of non-dimensional length of reinforcing beam to non-dimensional 
length of foundation beam is kept 1.5. The foundation beam is assumed to be smooth 
and no prestress forces are considered. 
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4.3.1 Ranges of the non dimensional parameters 

The ranges of the non-dimensional parameters used in the present 
study are tabulated in Table 4.2. 

Table 4.2: The ranges of non-dimensional parameters 


Non-dimensional parameters 

Ranges 

Non-dimensional depth of placement h 

0.5 -1.5 

Relative flexural rigidity of foundation beam and reinforcing 

beam, R = 

10-60 

k 

Relative stiffness of the soil layers, r = — 

1-20 

Non-dimensional unit weight of upper layer of soil,y = 

Qt 

0.05 - 1 


4.3.2 The effect of number of layers of reinforcement on Settlement response of 
foundation beam 

Following cases have been investigated for the study of effect of number of 
layer of reinforcement. 

a. Unreinforced case: This has been simulated using negligible bending stiffness 
of reinforcing beam i.e. E 2 I 2 is of the order 10' 

b. Single layer of reinforcement. 

c. Two layer of reinforcement. 

d. Three layer of reinforcement. 

e. Four layer of reinforcement. 

f. Five layer of reinforcement 
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All the above cases have been shown in Fig. 4.2 and Fig 4.3 for concentrated load 
at the centre of foundation beam and concentrated loads at the edges of foundation 
beam respectively. It is observed that increasing number of reinforcing layer beyond 
three does not have significant influence on the deflection profile of the foundation 
beam. The percentage reduction in the maximum normalized deflection of the 
foundation beam due to inclusion of various numbers of reinforcing layers with 
respect to unreinforced earth has been given in table 4.3 for load acting at centre and 
at edges respectively. 

Table 4.3: Effect of number of layers of reinforcements on maximum deflection 


of foundation beam 


Case 

Percentage reduction in maximum settlement with respect to 


unreinforced soil 


Concentrated load at centre 

Concentrated load at edges 

Single layer 

18.15 

22.76 

Two layer 

21.80 

30.34 

Three layer 

23.26 

34.83 

Four layer 

24.04 

37.93 

Five layer 

24.54 

40.00 


From Fig. 4.2 it is observed that for concentrated load at the centre of 
foundation beam, effect of number of layers of reinforcement on maximum 
normalized deflection i.e. at the centre of foundation beam is significant and the 
reduction in maximum normalized deflection is up to 23 % with respect to 
unreinforced earth when three layers of reinforcements are considered. It is found that 
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when one layer of reinforcement is used the percentage reduction in maximum 
normalized is maximum and is 18 % among all layers of reinforcement. The reduction 
in settlement due to adding second and third layer is 3-4 and 1-2 % respectively. 
Further adding reinforcement has negligible effect on deflection profile of foundation 
beam. 
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Fig. 4.2: Effect of number of layers of reinforcement on deflection profile of 
foundation beam subjected to concentrated force at its centre 
From Fig. 4.3 it is observed that the for concentrated load at edges of 
foundation beam, effect of number of layers of reinforcement on maximum 


normalized deflection i.e. at the edges of foundation beam is significant and the 
reduction in maximum normalized deflection is up to 35 % with respect to 
unreinforced earth when three layers of reinforcements are considered. This reduction 
is more than that of the case when load is acting at centre. The reason is that in case of 



load acting at edges foundation beam offers less resistance compared to the case when 
load is at the centre. Subsequently the reinforcing beam offers more resistance. It is 
found that when one layer of reinforcement is used the percentage reduction in 
maximum normalized deflection is maximum and is 23 % among all layers of 
reinforcement. The reduction in settlement due to adding second and third layer is 1 % 
and 5 % respectively. Further, adding more reinforcement has negligible effect on 
deflection profile of foundation beam. 



Normalized distance from centre (z) 

Fig. 4.3: Effect of number of layers of reinforcement on deflection profile of 
foundation beam subjected to concentrated force at its edges 
As it has been observed that increasing the number of layer of reinforcements 
beyond three does not have significant influence on the deflection profile of 
foundation beam, in the following sections for all the parametric studies only three 
layers of reinforcement are considered. For the study on deflection profile of 




reinforcing beams, the deflection profile for layer - 1 (the layer nearer to the 
foundation beam) and layer - 3 (the layer at the soil interface) are shown. The 
deflection profile for the middle reinforcing layer is not shown, and its deflection 
profile lies in between that of layer - 1 and layer - 3 . 

4.3.3 Effect of load 

The effect of load Q, is studied through varying non-dimensional 


parameter y 



)and keeping related parameters other than Q, constant. The 


parameter y decreases with increase in load Qt. Since Rj depends on parameters such 
as stiffness of foundation beam and coefficient of subgrade reaction of upper layer of 
soil layer, these are fixed as i? = 20, r = 5, Zr = 1.5, /?; = 1 and the same is used 
for further parametric studies in the following sections. 

In Fig 4.4 and Fig. 4.5, the deflection profiles for the foundation and the 
reinforcing beams are presented respectively for a set of typical values of parameters, 
showing the effect of a concentrated load at the center of the foundation beam through 
variation of y . Fig. 4.6 and Fig. 4.7 shows similar plots for the foundation beam 
subjected to end loads. The load parameter Q, appears in the non-dimensional 


deflection parameters y, = and ‘ 3 ' 

Q,R\ Qi^i 


in the denominator. Therefore, 


for lower loads higher non-dimensional deflections are observed which actually imply 
lower deflections in dimensional form. 

From Fig. 4.4 it observed that the maximum normalized deflection occurs at 
the centre and decreases to minimum towards the edge. The maximum normalized 
deflection at the centre decreases by 75 % and minimum normalized deflection at 
edge decreases from 0.253 to 5.32 as y decreases from 1 to 0.05. 
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Fig. 4.4: Normalized deflection of foundation beam for various values of y for 


load acting at the centre 

From Fig. 4.5 it is observed that the maximum normalized deflection occurs at 
the centre and decreases to minimum and becomes constant towards the edges. The 
maximum normalized deflections of layer — 1 and layer — 3 decrease by 76 % and 
77% respectively as y decreases from 1 to 0.05 and the minimum normalized 
deflections at the edge of layer - 1 and layer - 3 decrease by 97 % as y decreases 
from 1 to 0.05. It is seen that at lower load e.g. y = I layer - 1 has more influence 


than that of layer - 3. 
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Normalized distance from centre (z) 

Fig. 4.5: Normalized deflection of reinforcing beams for various values of y for 

load acting at the centre 

From Fig. 4.6 it observed that the maximum normalized deflection occurs at 
the edges and decreases to minimum towards the centre. The maximum normalized 
deflection at edge decreases by 75 % and minimum normalized deflection at the 
centre decreases by 84 % as y decreases from 1 to 0.05. 
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loads acting at the edges 

From Fig. 4.7 it is observed that the maximum normalized deflection occurs at 
z = 3.6 and decreases to minimum towards the edges. The maximum normalized 
deflections of layer - 1 and layer - 3 decrease by 78 % and 80 % respectively as y 
decreases from 1 to 0.05 and minimum normalized deflections at the edge of 
layer - 1 and layer - 3 decrease by 97 % as y decreases from 1 to 0.05. It is seen that 
at lower load y = 1 layer - 1 has more influence than that of layer - 3. The 
normalized deflections at the centre of layer - 1 and layer - 3 decrease by 85 % as y 


decreases from 1 to 0.05. 
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Fig. 4.7: Normalized deflection of reinforcing beams for various values of y for 


loads acting at the edges 
4.3.4 Effect of depth of placement of reinforcements 

In Fig 4.8 and Fig 4.9, the deflection profiles for the foundation and the 
reinforcing beams are presented respectively for a set of typical values of parameters, 
showing the effect of depth of placement of reinforcement for a concentrated load at 
the center of the foundation beam. Fig. 4.10 and Fig. 4.11 shows similar plots for the 
foundation beam subjected to end loads. 

From Fig. 4.8 it is seen that maximum deflection occurs at the centre of the 
beam and then decreases gradually to lowest towards the edges of the beam. With the 
increase in normalized depth of placement from 0.5 to 1.5, the maximum normalized 
deflection at the centre increases by 26 % while corresponding deflection at the edge 
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of the beam increases from 0.2 to 0.63. As the depth of placement increases all the 
reinforcement layers are offering less resistance resulting in increase in deflection. 



Fig. 4.8: Normalized deflection of foundation beam for various values of h’ for 


load acting at the centre 

Fig. 4.9, it can be seen that the maximum deflection of reinforcing beams 
occurs at the centre and decreases to a constant value towards the edge. The layer one 
deflects more than layer three for same depth of placement at the centre of 
reinforcement. This shows that reinforcing beam near to the foundation have greater 
influence. The increase in maximum deflection for layer - 1 and layer - 3 at centre of 
beam is around 27 % and 30 % respectively as normalized depth of placement 
increases from 0.5 to 1.5. The minimum normalized deflection at the edge increases 
by 80 % for both layer -1 and layer - 3 as h increases from 0.5 to 1.5. Deflection 
profiles for various normalized deflection become parallel at the edges. This shows 
that deflections at edges are mainly due to surcharge. 
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Fig. 4.9: Normalized deflection of foundation beam for various values of h’ for 


load acting at the centre 

From Fig. 4.10 it is observed that maximum deflection in the foundation beam 
occurs at the point of application of the load, i.e., at the edge of the beam, then 
decreases gradually towards the centre of the beam and becomes least at the centre of 
the beam. It is observed from the figure that as the normalized depth of placement 
increases from 0.5 to 1.5, the maximum normalized deflection at edge of foundation 
beam increases by 25 % and normalized deflection at the centre increases by 48 %. 
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Fig. 4.10: Normalized deflection of foundation beam for various values of h’ for 


loads acting at the edges 

Fig. 4.11 shows that the maximum deflection occurs at z = 3.6, decreases 
rapidly towards the edge of beam because of less resistance from foundation beam 
and gradually decreases towards the centre of beam and becomes constant. The 
minimum deflection occurs at the edges of reinforcing beams. As normalized depth 
of placement increases from 0.5 to 1.5 the maximum normalized deflections at z = 3.6 
increase by 32 % and 35 % for layer - 1 and layer - 3 respectively. The normalized 
deflections at the centre increase by 50 % and 53 % for layer - 1 and layer -3 
respectively. 
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Fig. 4.11: Normalized deflection of reinforcing beams for various values of h’ for 


loads acting at the edges 
4.3.5 Effect relative flexural rigidity of beams 

In Fig 4.12 and 4.13, the deflection profiles for the foundation and the 
reinforcing beams are presented respectively for a set of typical values of parameters, 
showing the effect of relative flexural rigidity of beams for a concentrated load at the 
center of the foundation beam. Figs. 4.14 and 4.15 show similar plots for the 
foundation beam subjected to end loads. The relative flexural rigidity R, affects the 
settlement response of foundation beam and reinforcing beams significantly. 
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From Fig. 4.12 it is observed that for the foundation beam the maximum 
deflection occurs at the centre of the beam and gradually decreases to minimum at the 
edge. The straight portion of the deflection profile towards the edge shows that the 
deflection at the edge is mainly due to rigidity of the beam. It is observed that with the 
increase in the ratio R there is an increase in maximum deflection and this increase in 
deflection becomes negligible after a value of R = 40. As the value of R increases 
from 10 to 60 the maximum normalized deflection increases by 4 %. For higher 
values of R (R = 40 and R =6 0) negative deflections at the edges are observed. 



Fig. 4.12: Normalized deflection of foundation beam for various values of JR for 

load acting at the centre 

Fig.4.13 shows that the maximum normalized deflection is observed at the 
centre of the beams, it gradually decreases towards the edge and finally becomes 
constant. It is observed that as ratio R increases from 10 to 60 the normalized 
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deflections decrease by 82 % and 66 % for layer - 1 and layer - 3 respectively. As 
ratio R increases flexural rigidity of reinforcing beam decreases hence higher 
deflections are observed. The minimum normalized deflections at the edges decrease 
by 73 % for both layer - 1 and layer -3. 



Fig. 4.13: Normalized deflection of reinforcing beams for various values of R for 

load acting at the centre 

Fig. 4.14 shows that the maximum deflection is observed at the edges of the 
beam and it increases by 17 % as the ratio R increases from 10 to 60. The minimum 
deflection is observed at the centre of the beam and it decreases by 7 % as the ratio R 


increases from 10 to 60. 
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Fig. 4.14: Normalized deflection of foundation beam for various values of R for 

loads acting at the edges 

From Fig. 4.15, it is seen that the maximum normalized deflection occurs at 
z = 3.6. It is observed that as the ratio R increases from 10 to 60 maximum normalized 
deflections at z =3.6 decrease by 80 % for both layer - 1 and layer -3. The normalized 
deflections at the centre decrease by 85 % for both layer - 1 and layer - 3. The 
minimum normalized deflections the edges decrease by 96 % for both layer - 1 and 


layer -3. 
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Fig. 4.15: Normalized deflection of reinforcing beams for various values of i? for 

loads acting at the edges 

4.3.6 Effect of relative stiffness ratio of the soils 

In Fig 4.16 and Fig. 4.17, the deflection profiles for the foundation and the 
reinforcing beams are presented respectively for a set of typical values of parameters, 
showing the effect of variation of stiffness ratio of upper layer soil and lower layer 
soil r for a concentrated load at the center of the foundation beam. Figs. 4.18 and 
Fig. 4.19 show similar plots for the foundation beam subjected to end loads. 

From Fig. 4.16 it is observed that for the load acting at the center, the 
maximum normalized deflection on foundation beam occurs at the centre and 
increases from 0.7 to 5.24 as the ratio r increases firom 1 to 20. As the ratio r increases 
the lower soil becomes poor and hence the increase in deflection is observed. The 
minimum deflection at edge increases from 0 to 2.57 as r increases from 1 to 20. 
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Fig. 4.16: Normalized deflection of foundation beam for various values of r for 

load acting at the centre 

Fig. 4.17 shows that the maximum nonnalized deflection, occurs at the centre 
of the beam, increases from 0.03 to 0.26 and 0.018 to 0.25 for layer -1 and layer - 3 
respectively as r increases from 1 to 20 and gradually decrease to minimum deflection 
at the edge of the beam. The normalized deflection at edge increases from 0.004 to 
0.023 as r increases from 1 to 20 for both layer - 1 and layer -3. 
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Fig. 4.17: Normalized deflection of foundation beam for various values of r for 


load acting at the centre 

Fig. 4.18 shows that the maximum normalized deflections are found at the 
edges and it increase from 0.871 to 5.1 and minimum normalized deflection at the 
edge increases from 0.33 to 4.4 as r increases 1 to 20. 
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Fig. 4.18; Normalized deflection of foundation beam for various values of r for 


loads acting at the edges 

From Fig. 4.19 it can be seen that the maximum deflection occurs at z = 3.6 of 
foundation beam and minimum settlement at the edge of reinforcing beam. It is also 
observed that for a typical value of ratio r layer one deflects more than the layer three. 
It shows that beam near to the foundation has greater effect. The maximum deflection 
increases 0.028 to 0.239 and 0.0163 to 0.229 for layer - 1 and layer -3 respectively as 


r increases from 1 to 20. 
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Fig. 4.19: Normalized deflection of foundation beam for various values of r for 


loads acting at the edges 


4.4 Particular case of multiple layer reinforcement without 
reinforcement at interface 

As a particular case of the above study, only one reinforcement layer is kept at 
the mid depth of dense compacted soil layer. This is simulated by analyzing the two 
layer reinforcement of the present model and keeping very low value of E 2 I 2 very low 
i.e. of order 10'^° for the beam elements corresponding to the reinforcement at the 
interface of soil layers. 

In Fig 4.20 and Fig 4.21, the deflection profiles for the foundation and the 
reinforcing beams are presented respectively for a set of typical values of parameters, 
showing the effect of reinforcement at the mid depth of dense compacted sand layer 
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for concentrated load acting at the center of the foundation beam. Figs. 4.22 and 
Fig. 4.23 show similar plots for the foundation beam subjected to end loads. 

It is observed that the reduction in the maximum deflection of foundation 
beam is 20 % and 30 % with respect to unreinforced for load acting at centre and edge 
respectively. 



Fig. 4.20: Normalized Deflection of foundation for reinforcement at the mid 
depth for load acting at the centre 
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Fig. 4.21: Normalized Deflection of reinforcement for reinforcement at the mid 
depth for load acting at the centre 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Normalized distance from centre (z) 

Fig. 4.22: Normalized deflection of foundation for reinforcement at the mid 
depth for load acting at the edges 



iMormaiizea aeTiection of reinforcing beams 


R = 20, r = 5, R, = 1 , ki = 100MN/m3, ^=0.08, h'=1 .Q* = 100000N 



Fig. 4.23; Normalized deflection of reinforcing beam for reinforcement at the 
mid depth for load acting at the edges 



CHAPTERS 


CONCLUSIONS AND SCOPE OF FUTURE STUDY 

5.1 General 

Following the problem formulation and method of analysis, multilayered 
reinforced earth foundation system is analyzed for the settlement response of 
foundation and reinforcing beams. The present model results are validated and then 
detailed results for parametric studies are obtained and discussed. Based on the results 
of the present model study following conclusions are drawn. 

5.2 Conclusions 

1. Results of the present model are compared with the close form solution by 
Maheshwari (2004) for a degenerated case. The present model results are in 
very good agreement with close form solution by Maheshwari (2004). 

2. The number of layers of reinforcement significantly affects the settlement 
response of foundation beams for load acting at centre as well as at edges. The 
reduction in maximum normalized deflection for load acting at centre and for 
loads acting at the edges are 23 % and 35 % when three layer of 
reinforcements are used. Number of layers of reinforcement more than three 
does not have any significant effect on the settlement response of foundation 
and reinforcing beams for both the loading conditions. 

3. The Load acting on the foundation beam affects the settlement response of 
both foundation and reinforcing beams significantly. For the load acting at 
center, the decrease in the maximum normalized deflection at the centre for 
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foundation beam and reinforcing beams is of the order of 75 % and 76 - 77 % 
when Y decreases from 1 to 0.05. For loads acting at edges, the decrease in the 
maximum normalized deflection of the foundation beam and reinforcing 
beams are 75 % and 78 - 80 % respectively as y decreases from 1 to 0-05. 

4. Normalized depth of placement affects the settlement response of foundation 
beam and reinforcing beams significantly for the both the loading conditions. 
For the load acting at centre the increase in maximum deflection at the centre 
for foundation beam and reinforcing beams is 26 % and 27 - 30 % when 
normalized depth of placement increases from 0.5 to 1.5. For the loads acting 
at the edges the maximum normalized deflection of foundation beam and 
reinforcing beams increases by 25 % and 32- 35 % respectively when 
normalized depth of placement increases from 0.5 to 1.5. 

5. The relative flexural rigidity ratio R affects the settlement response of 
foundation beam and reinforcing beams significantly for both the loadings. 
For the load acting at centre the maximum normalized deflection at the centre 
of foundation beam increases by 4 % when R increases from 10 to 60. The 
maximum normalized deflection of reinforcing beams decreases by 60 - 90 % 
when R increases from 10 to 60. For the loads acting at edges the increase in 
maximum deflection at the edge of foundation beam is 17 % when R increases 
from 10 to 60 and the maximum normalized deflection of reinforcing beams 
decreases by 80 % when R increases from 10 to 60. 

6. The relative stiffness ratio r of the soils affects the settlement response of 
foundation beam and reinforcing beams significantly for both the loading 
conditions. For the load acting at centre the increase in maximum normalized 
deflection at the centre of foundation beam and reinforcing beams are 650 % 
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and 700 % respectively when r increases from 1 to 20. For the loads acting at 
edges the maximum normalized deflection of foundation beam increases by 485 
% when r increases from 1 to 20 respectively. 

7. When only one layer of reinforcement is placed at the middle of compacted soil 
layer the reduction in maximum normalized deflection with respect to 
unreinforced earth is 20 % and 30 % for load acting at centre and loads acting at 
the edges respectively. 

5.3 Scope of future study 

Following are the future scopes of the present work 

1. Higher order approximation for foundation beam and reinforcing beam are to 

be taken, — ^ and — ^ are to additional degrees of freedom so that shear 
dx^ dx^ 

force distribution and bending moment distribution on beams can be 
determined. 

2. The present model results are to be validated by experiments. 

3. Dynamic consideration is to be given to consider the inertial effect and study 


of moving loads. 
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